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Abstract 



We have determined the phase diagram of the simplest version of a lattice model 
introduced in the recent work of Kazakov and Migdal. If mo and A are the bare 
mass and self coupling of the scalar field in the model respectively, we find a line of 
first order phase transitions in the (mo,A) plane ending in a critical point where A is 
nonzero. Kazakov and Migdal speculate that their model of scalar field theory could 
induce QCD. Our work indicates that for N = 2 there is no continuum limit for the 
Kazakov-Migdal model except at the critical end point. Whether or not a nontrivial 
continuum limit exists in the vicinity of the critical point requires careful study of the 
renormalization group properties of the model. 



Recently Kazakov and Migdal (KM) []T| came up with an intriguing idea: consider the 
following model, defined on a d-dimensional hypercubic lattice with action 



s=ny: 



(1) 



where $(a;) is a traceless scalar field in the adjoint representation of SU{N), covariantly 
coupled to the gauge fields U^{x) defined on the links of the lattice. Kazakov and Migdal 
speculate that perhaps, when the parameters of the scalar potential V are tuned appropri- 
ately, the model Eq. (1) will have a non-trivial continuum limit. Moreover the resulting 
continuum theory might actually be equivalent to QCD. 

If the above were true this apparent reformulation of QCD should not be considered as 
merely an interesting curiosity. The potential of this speculation lies in the fact that Eq. (1) 
in many ways is much simpler than conventional lattice QCD. As a matter of fact, based on 
the work of KM and of Migdal it looks as though in the limit of large the model might 
be analytically tractable and that the degrees of freedom associated with the eigenvalues of 
$ could actually be the long sought "master field" of large A^ QCD. 

In order to see how such a scenario could possibly be true and also to motivate the 
present study, let us examine Eq. (1) in more detail. Integrating out the scalar fields in the 
partition function will give rise to an induced action for the gauge fields which because of 
the gauge invariance of Eq. (1), will be given as an (infinite) sum over Wilson loops on the 
lattice. Moreover, because the scalar field is in the adjoint representation, the Wilson loops 
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will also be. Concretely, for V{4>) = ^trcf?, the integral is Gaussian and one obtains the 
induced action 

2'(r)-i|tr[/(r)|2 
r /[r]mo 

where 1{T) is the length of the path F. Note that Eq. (2) inherits from Eq. (1) a local center 
symmetry f/^(x) — ^ Z^{x)U^{x), a fact recently discussed by Kogan, Semenoff and Weiss [Q. 
We will return to this point again later on. 

How could the theory defined by Eq. (2) possibly be related to QCD? To begin with 
the induced action in Eq. (2) is as good a lattice action as any other in the sense that in 
the naive continuum limit, it produces a term proportional to F^^^. If a continuum limit 
indeed exists for this model, we expect that only this term will survive based on the familiar 
renormalization group argument. The model is unusual only because the basic building 
blocks of the action are arbitrary loops ( in the adjoint representation). In the weak smooth 
gauge field limit this model has a critical point at ttIq = 2d as can be seen easily from 
Eq. (1) where the lowest order action becomes an action for a free scalar field. Taking 
fiuctuations into account the critical point is expected to shift to a different value m^. Using 
the definition = mg — KM calculated the coefficient of the F^^ term and obtained 



Omd — 2^ j^^^ 2Z{r) ' 



the correspondence 
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in the smooth gauge fields hmit. 



KM and Migdal [|| pointed out that the model in Eq. (1) can be easier dealt with 
analytically and numerically if one first integrates out the gauge fields. In the large limit, 
the resulting effective action for the master field is best rewritten as a two-matrix model, 
which can be treated using modern matrix model technology. KM suggest that the master 
field gives nontrivial scaling laws for physical quantities (for example, glueball mass Mg) 



Mi 



(m 



2n7 
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where the value of 7 is different from the usual trivial scalar theory (7 = 1 with possible 
logarithmic corrections). KM have given the following argument for the value of 7: if the 
continuum theory of Eq. (1) is to "induce" QCD correctly, we must view the scalar as a 
heavy "constituent" field whose mass acts as an effective UV cutoff for the gauge field in 
the continuum. This immediately gives the usual relation between a physical quantity , the 
cutoff and the bare coupling 

(5) 
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Combining this expression with that of Eq. (|^) one gets the scaling relation for Mg and find 
7 = 23/22 (This is the value in the smooth gauge field limit. When the gauge field is treated 
nonperturbatively 7 value can be quite larger than one 0). 

Based on the above discussions we expect that in order for Eq. (1) to correctly induce 
QCD, we need two things, a critical point and nontrivial scaling of physical quantities. In 
the present letter we follow KM's suggestion and study the simplest version of the model, 
SU{2) (A^ = 1 is trivial). Certainly the possibility of Eq. (1) inducing QCD does not 
crucially depend on being large. Integrating out the gauge fields one obtains the following 
partition function |I| 



<t>>o 



[dcj)] exp 




2V(0) 



+ ^ /n 

<xy> 



sinh{4:(f){x)(f){y)) 
(j){x)(j){y) 
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where < xy > denotes the sum over nearest neighbors of x. The field 0(x) in Eq. (^ is the 
(positive) eigenvalue of the matrix $ and the restriction to (f){x) > expresses the fact that it 
is a "radial" variable. According to what was said before, we must find a critical point close 
to which a continuum limit can be constructed. It is well known that the continuum limit 
of a scalar model with polynomial interaction terms is trivial with vanishing renormalized 
coupling . In principle the continuum limit of Eq. (P) could avoid the triviality problem 
due to the presence of the nonpolynomial interaction terms (if we could integrate out the 
angular variables in an 0{N) model exactly, however, the resulting effective action for the 
radial variables would be also nonpolynomial). Whether or not this is the case and Eq. (^ 
allows for nontrivial scaling is numerically a much more difficult question and remains to be 
answered in the future. 



We restrict ourselves to a scalar potential of the form V{(f)) = ^trcj)^ + 



A nonzero 



A is necessary to stablize the system and allows numerical simulation in the small itlq region. 
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In order to get an idea of what the phase diagram looks hke in the (A, mo) plane, we first use 
a simple mean field theory (MFT), i.e. we set = and assume u{x) = constant. We 
then analyzed the resulting potential looking for minima and phase transitions. As it turns 
out the qualitative features of the phase diagram from the MFT effective potential analysis 
are in excellent agreement with Monte Carlo simulations. In Fig. 1, we show a typical plot 
of the MFT potential Veff{u) at a rather small value of A = 0.01. There obviously is a first 
order phase transition at mg = ml = 7.65, which is not very far away from the smooth gauge 
field limit ml = 8. For this particular value of A we show a comparison of MFT and Monte 
Carlo for 0-field vacuum expectation value in Fig. 2. All our Monte Carlo simulations were 
done on 8^ lattice using a single-hit Metropolis algorithm on the field u. The errors in Fig. 
2 were obtained by blocking averages and < > values are the result of 5K sweeps through 
the lattice at each point. Clearly, the agreement between MFT and Monte Carlo is very 
good. The transition is obviously first order with a very clear hysteresis in < > (we also 
have obtained clear two state signals at m^). Interestingly then, our first conclusion is that 
in the A — limit no continuum limit exists. 

The appearance of first order phase transitions does not come as a surprise. The simple 
SU{2) lattice theory in the adjoint representation {SU{2)/Z2 = 5*0(3)) has action S = 
PaY^d \TrUj:,\'^ and possesses a first order phase transition at some critical value of (3a 0- 
Presumably this transition is related to the Z2 invariance of the action 0. The question now 
is whether a first order transition can be avoided by tuning A. The answer is yes: In MFT 
the gap in -u (or < >) decreases rapidly as A is increased. For A > 2.57 the gap disappears 
all together and there is no phase transition anymore. Thus the phases of strong and weak 
mo are analytically connected. The point (A = 2.57, mg = 4.515) is a critical point in MFT: 
V^'j-j(0) vanishes at the minimum of the potential. In Fig. 3 we show how this MFT scenario 
is confirmed by Monte Carlo simulations. The location of the critical endpoint by Monte 
Carlo is not very precise. It is taken at the value of A where the hysteresis effect in < > 
disappears and the critical value of m^ is determined by the position of the peak of (/)-field 
susceptibility 

X = TT E [< 'P{^)'P{y) > - < <P{^) >< <P{y) >] • (7) 

x,y 

In conclusion, we have found a line of first order phase transitions in the (A, mp) plane 
of the SU(2) KM model. There is no continuum limit at A ^ 0. The first order phase 
transition line terminates in a critical point at A > 0, in the vicinity of which one may 
attempt to construct a continuum theory. The precise nature of this continuum theory 
remains to be determined. This is a challenge for further research. 
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Figure Caption 



Figure 1: The MFT effective potential Veff{u) at A = 0.01. The two minima become 
degenerate at ml = 7.65 which agrees with the observed critical point for the first order 
phase transition in Monte Carlo simulation. 

Figure 2: Comparison of the mean field theory and Monte Carlo simulation result for 0-field 
vacuum expectation value. The solid hues are the MFT predictions and Monte Carlo data 
points are indicated by *. The errors are smaller than the size of the symbols. 

Figure 3: The phase diagram. The solid line is the mean field theory prediction and 
the simulation results are indicated by diamonds. The end point of the first order phase 
transition line is marked by * for MFT and square for Monte Carlo estimate. At the end 
point the theory becomes critical. 
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